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Abstract. Temperature contrasts and magnetic field strengths of sunspot umbrae broadly fol- 
low the thermal-magnetic relationship obtained from magnetohydrostatic equilibrium. Using a 
' compilation of recent observations, especially in molecular bands, of temperature contrasts of 
starspots in cool stars, and a grid of Kurucz stellar model atmospheres constructed to cover 
O ' layers of sub-surface convection zone, we examine how the above relationship scales with ef- 
fective temperature (T e //), surface gravity g and the associated changes in opacity of stellar 
photospheric gas. We calculate expected field strengths in starpots and find that a given rela- 
tive reduction in temperatures (or the same darkness contrasts) yield increasing field strengths 
against decreasing T e ff due to a combination of pressure and opacity variations against T e //. 
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2 ■ 1- Introduction 

Despite a lack of deductive magnetohydrodynamic explanation for the formation and 
the equilibrium of a sunspot, extensive observations in combination with magnetohy- 
drostatic models have provided reasonable understanding of the thermal-magnetic struc- 
ture of sunspots in the observable layers. Reduced temperature and gas pressure in- 
side sunspots dictate dominantly the thermal-magnetic relationship derived from the 
magneto-hydrostatic balance. Such a relationship causes the Wilson depression - geo- 
metrical depression of the observable optical depth unity level within a sunspot -, which 
provides an explanation for the field strengths observed in sunspots and also relates the 
intensity (or the brightness) contrasts to field strengths. Here we examine how such a 
thermal-magnetic relationship scales with the stellar parameters, viz. the effective tem- 
perature T e f f and surface gravity g as well as the associated changes in the opacity of the 
stellar photospheric gas. We then discuss the implications of such scalings for the inter- 
pretations of observed field strengths. We also discuss how such scalings could be crucial 
players in the activity related photospheric brightness variations and their correlation 
with other activity measures. 



2. Method of Calculation 

Thermal-magnetic relationship for starspots is obtained from a simplified magnetohy- 
drostatic (MHS) condition, that relates the magnetic field strength to the temperature 
at the axis of a vertical column of magnetic field (Solanki et al. (1993)). The radial com- 
ponent of the magneto-hydrostatic force-balance equation (Maltby 1977, Solanki et al. 
1993), after neglecting the magnetic curvature force and with the use of the equation of 
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state P = RpT / /i, yields the thermal-magnetic relation, 



T(r, z) = p(r,z)p e (z) 
T e (z) n e (z)p(r,z) 
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(2.1) 



where P, T, p, p and R denote the gas pressure, temperature, density, mean molecular 
weight and the gas constant respectively. The subscript e stands for the external atmo- 
sphere. B z is the vertical component of magnetic field and r is the radial distance from the 
center of the spot and z is the depth measured from continuum optical depth r Cje =l. We 
further neglect the r-dependence of quantities in the above equation, thus the calculated 
quantities refer to the axis of the spot, and we refer B z simply as B hereafter. Eqn l2.ll 
is valid for each level z. The variation of external atmospheric quantities with depth z 
are determined by g and T e / / of the parent star and are taken from the grid of Kurucz 
stellar models constructed to cover deeper regions of the convection zone using the AT- 
LAS9 stellar atmosphere code (Kurucz 2001). We prescribe temperature contrasts for 
starspots under two cases, case (i): use the empirical relation AT = (590. * logg) — 680!^ 
(O'Neal et al. (1996)1 to determine the effective temperature of the spot T e ff ySpt in a 
parent star characterised by g and T e ff. T spt (T c = 2/3) = T e ff >sp t = T e ff — AT, case 
(ii): the temperature contrasts are independent of g and are of a constant ratio of T e ff. 
AT = 0.3T e // (Berdyugina (2005)). Hydrostatic equilibrium inside the spot yields, to a 
very good approximation (Cox and Giuli 1968), 



P(t c = 2/3) 
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The density inside the spot at r c = 2/3 is determined by solving the above equation 
with the use of Rosseland mean opacities kr from the tables of Kurucz (1993) and 
Alexander and Ferguson (1994). Saha's equation is used to determine the mean molecular 
weight. The density difference between the r c = 2/3 and r c = 1 levels within the spot 
is assumed to be negligible, i.e., p(r c = 1) w p{r c = 2/3). Using the assumption that 
p(r c = 1) = p e (z — Z w ), i.e., the densities inside and outside the spot are equal at the 
level of Wilson depression Z w (Maltby (1977)), we find the depth z at which the above 
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Figure 1. Variation of field strength B (solid lines) and the Wilson depression Z w (dashed 
lines) against log(g) for T e // = 6000A' 
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equality is satisfied in the parent stellar model. This gives the value of Z w , and now we 
determine the only unknown quantity B at this level using Eqn l2.ll 

3. Results and Discussions 

Figures 1 and 2 summarise the main results of our study. Each figure shows the vari- 
ation of observable field strengths B(Z W ) on the spot axis as a function of log(g) for 
a particular T e // for the two distinct kinds of variation of spot temperature given as 
case (i) (red curves) and (ii) (green curves) in the previous Section. The Wilson depres- 
sions, Z w , are shown as dashed curves and the B values as solid curves. Z w are scaled in 
terms of pressure scale heights H p at the r c =l level in the quiet photospheres. Fig.l, for 
T e ff = 6000-ftT, shows that the B values are not very different for the two cases above, 
and moreover the Wilson depressions Z w are of the order of scale heights. In contrast, 
Figure 2, for T e f j = 5000X , shows that the results for the above two cases are very differ- 
ent: the gravity dependent temperature contrasts (as implied by the observed relation of 
O'Neal et al.(1996)) requires very strong magnetic fields for starspot equilibrium, which 
stems from the large values of Z w . On the other hand, a constant value of temperature 
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Figure 2. Same as Fig. 1 but for T e ff = 5500A" (top) and 5000^ (bottom) 





122 Rajaguru and Hasan 

reduction (case ii above) over the range of surface gravity values used yields B values 
within the observed ranges (see Berdyugina (2005) and references therein). Consequently, 
a relatively less temperature reduction and hence a less amount of gas evacuation is suf- 
ficient to attain a given value of field strength in cooler stars than that in stars hotter 
than about 5500/C In other words, spots of a given field strength would appear much 
darker in a star of T e tf = 6000if than those in T e ff = 5000iC 

According to Lockwood et al.(1992), younger and faster rotating stars are 'spot domi- 
nated', i.e. more flux in spots than small-scale fields and faculae, and hence grow darker 
as activity increases. If such a 'spot domination' is purely dependent on age (rotation) 
but not on color (spectral type), then it could be thought of as a phenomenon not con- 
tradictory with the g dependent temperature contrasts for spots derived by O'Neal et al. 
(1996). However, our results in Figure 2 imply unrealistically large B and Z w values for 
spots in such cases. On the other hand, consistency between Lockwood et al.'s results 
requiring 'spot domination' for younger stars and the situation of case (ii) of our results 
requires that there be a color T e f f dependence of spot properties in addition to the age 
dependence. 

Alternatively, results of Lockwood et al.(1992), but without the requirement of 'spot 
domination' in younger stars, could be made consistent with our case (ii) of less gas 
evacuated spots if the small-scale fields forming the 'faculae' too are of such less evacu- 
ated state. This would imply that faculae in younger and cooler stars are less bright and 
therefore these stars grow darker as activity increases. Interestingly, the superadiabatic- 
ity that drives the convective collapse of small-scale flux tubes indeed linearly decrease 
with T e f f and is found not to intensify such fields to a high degree of evacuation as in 



stars hotter than about 5000 K (Rajaguru et al. (2002)). Hence, it would appear that 



the near-surface thermal structure of cool stars crucially determine the key properties of 
magnetic structures small and large. We conclude that the scaling of thermal and mag- 
netic properties of starspots with both the gravity g (age) and T e ff (color) are crucial 
for a consistent interpretation of observed correlations between activity measures that 
sample the different heights in the outer atmospheres. 
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